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We calculate the magneto-optical conductivity and electronic density of states for silicene, the 
silicon equivalent of graphene, and similar crystals such as germanene. In the presence of a perpen- 
dicular magnetic field and electric field gating, we note that four spin- and valley-polarized levels 
can be seen in the density of states and transitions between these levels lead to similarly polarized 
absorption lines in the longitudinal, transverse Hall, and circularly polarized dynamic conductivity. 
While previous spin and valley-polarization predicted for the conductivity is only present in the 
response to circularly polarized light, we show that distinct spin- and valley-polarization can also 
be seen in the longitudinal magneto-optical conductivity at experimentally attainable energies. The 
frequency of the absorption lines may be tuned by the electric and magnetic field to onset in a range 
varying from THz to the infrared. This potential to isolate charge carriers of definite spin and valley 
label may make silicene a promising candidate for spin- and valley tronic devices. 

PACS numbers: 78.67.Wj, 78.20.Ls, 71.70.Di, 72.80.Vp 



Introduction: Silicene, a monolayer of silicon atoms 
bonded together on a two-dimensional (2D) honeycomb 
lattice, has recently been synthesized^ and begun to 
garner increased theoretical attention 7 - - — as it features a 
Dirac-like electron dispersion at the K points of the Bril- 
louin zone and promises to exhibit exciting properties 
beyond those present in graphene. The similarities with 
graphene result from carbon and silicon residing in the 
same column on the chemical periodic table. The larger 
ionic size of silicon atoms causes the 2D lattice of silicene 
to be buckled^ ' 10 ' 12 such that sites on the A and B sub- 
lattices sit in different vertical planes with a separation 
of d ~ 0.46 A 11 ' 12 as illustrated in Fig. [TJ Because of the 




FIG. 1. (Color online) The crystal structure^ of silicene is 
based on the honeycomb lattice; the ionic size of the silicon 
atoms causes the A (red) and B (blue) sublattices to sit in 
separated parallel planes. 



buckled lattice, an onsite potential difference (A z — E z d) 
arises between the A and B sublattices when an electric 
field E z is applied perpendicular to the plane. Silicene is 
also predicted to have a stronger intrinsic spin-orbit gap 
than is seen in graphene^ with values (which can be in- 
creased under strain2ii&) predicted to be A so » 1.55 meV 



by density functional theory calculations ^ 10 ' 12 and is 
quoted as A so « 7.9 meV in tight-binding calculations 10 . 
The larger spin-orbit interaction makes silicene suscepti- 
ble to spin manipulation. The resulting band gap near 
the two valleys K and K'(= —K) of the first Brillouin 
zone provides a "mass" to the Dirac electrons that can be 
controlled by the strength of A ^ 11 ' 12 ' 14 . Silicene has also 
been predicte d 12 ' 14 to undergo a transition from a topo- 
logical insulator (TI) (an insulator with a gapless spec- 
trum of edge state s 21 ' 22 ) to a band insulator (BI) as the 
strength of A 2 becomes greater than A so . These qualities 
are also predicted for a monolayer of germanene which is 
isostructural to silicene but is expected to exhibit a much 
stronger spin-orbit band gap of A so w 24 — 93 meV from 
first principles^ and tight-binding calculations^. 

When subjected to an external magnetic field, Lan- 
dau levels (LLs) form in the electronic density of states 
and transitions between these levels result in absorption 
lines in the optical conductivity. This has been discussed 
theoretically^— and confirmed experimentally 2 ^— for 
graphene. In graphene, the n = LL is pinned at e = 
q ^-25,34-37 when an excitonic gap is included^^, the 
n = level splits between two distinct valley-polarized 
spin-degenerate energies but spin-polarized charge car- 
riers are not expected in this case. Conversely, due to 
spin-orbit coupling (SOC) and the response of silicene to 
an external perpendicular electric field, spin- and valley- 
polarized charge carriers appear due to the n — LL 
splitting between four spin and valley distinct energies. 

In this Letter, we examine the effect of exposing sil- 
icene (or similar 2D crystals, such as germanene 1 ^) to 
both external magnetic and electric fields with partic- 
ular attention to the valley- and spin-polarized regions 
of the electronic density of states and dynamical con- 
ductivity. This is of particular technological interest as 
four distinct valley- and spin-polarized currents can be 
generated through use of an in-plane electric field to gen- 
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erate a Hall current of spin- or valley-selected charges on 
the edges of the sampl o 38,39 . This ability to control the 
spin and valley index for use in such applications as data 
storage and data transmission is integral to valle y 38 i 40 ~ ^2 
and spintroni o 43 ' 44 devices. 

Low Energy Hamiltonian: It has been shown that the 
low-energy physics of silicene can be captured by the sim- 
ple low-energy Hamiltonia n 10 ' 13 ~ — 
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where £ = ±1 at the K and K' points, respectively, a z 
is the Pauli matrix associated with the electron's spin, 
t are the Pauli matrices associated with the sublatticc 
pseudospin, v w 5 X 10 5 m/s is the Fermi velocity of sil- 
icene and p x and p y are components of the momentum 
measured relative to the K/K' point. The first term 
of the Hamiltonian is the familiar graphene-type Dirac 
Hamiltonia n 45 ' 46 . The second term is the Kane-Mele 
term for intrinsic SOG^i while the final term is associated 
with the onsite potential difference between the two sub- 
lattices that results from an external electric field 1 - 3 - - — . 
There is also a Rashba SOC but it may be ignored 1 ^ 
for our purpose as it is a factor of 10 smaller than the 
intrinsic SOC. 

If we consider a perpendicular magnetic field and work 
in the Landau gauge, the magnetic vector potential, A = 
V x B, is written as A = (—By, 0,0). We make the 
Peierls substitution pi —> pi + (e/c)Aj in Eq. |T]) and 
obtain the low energy dispersion 



e (A ) = J sgn(n)^Aj a + 2\n\v 2 heB/c, n = ±l,±2,.. 

(2) 

where A^ a = — icr£A so + ^A 2 , a = ±1 for spin up 
and down, respectively, e is the elementary charge and 
B is the strength of the external magnetic field. The 
Zeeman energy is small and can be ignored 1 ^ 15 ' 47 . In 
Fig. [5] we show a schematic of the spin split bands at 
the K and K' point for B = with the LLs for finite 
B shown as dots. The form of the zeroth LL has strong 
implications when considering the LLs in the TI (A z < 
A so ) vs. BI (A z > A so ) regimes which can be seen in 
the upper frame (TI) and lower frame (BI) of Fig. [21 
respectively. Given our expression for £o(A^ CT ), we can 
see that at both K points in the TI regime, the n = 
spin up LL is at positive energy while the n = spin 
down level is at negative energy; refer to the upper frame 
of Fig. [2] In the BI regime, A so is now greater than A z 
and, therefore, the signs of Eo(^-Kf) an< ^ £ q(^-K'i) switch 
so that the n = spin up level is below zero at the K 
point while the spin down n = level is above zero at 
the K' point as illustrated in the lower frame of Fig. [5] 
This LL shift arises from the band inversion that results 
from the transition between the TI and BI regimes^. 

Electronic Density of States: Plots of the single-spin 
electronic density of states in the TI and BI regimes are 
shown in Fig. [31 There are four spin-polarized peaks 




FIG. 2. (Color online) Schematic of the LLs (solid dots) com- 
pared to the band structure for B = (solid curves). Blue/red 
represents spin up/down. The numbers identify the LL index. 
Upper: TI regime (A« < A so ) at the K (left) and K' (right) 
points of the Brillouin zone. Lower: BI regime (A z > A so ). 
A chemical potential /i (dashed green line) is shown that al- 
lows for spin- valley-polarized absorption lines in the magneto- 
optical conductivity due to the transitions identified with ar- 



which onset at to = —Aki, Ak>i, —Ak^ and Ak>-\ cor- 
responding to the four n — LLs. The remaining higher 
energy peaks (not shown) are spin degenerate. As scan- 
ning tunnelling spectroscopy experiments have been suc- 
cessful in observing LLs in graphene 4 ^— , similar work 
on silicene should detect four distinct low energy spin- 
and valley-polarized levels. 

Dynamical Conductivity: Using the standard Kubo 
formul a 16 ' 51 " —, we can derive expressions for the longitu- 
dinal and transverse Hall dynamical conductivities which 
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also yield the familiar selection rule; 
for LL transitions. For the absorptive part of the longi- 
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FIG. 3. (Color online) Spin-dependent density of states (red: 
spin down and blue: spin up). Left: in the TI regime 
(A s = 0.5A SO ). Right: in the BI regime (A z = 1.5A so ). 
Each plot is shown for B/A 2 o =0.0657T/meV 2 with a broad- 
ening parameter of r\ — 0.07A so . There are four spin- and 
valley-polarized peaks at low energy with all higher peaks 
(not shown) being spin degenerate. 



tudinal conductivity we find 
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where do = e 2 /(Ah) is the universal background conduc- 
tivity of graphene. For the absorptive imaginary part of 
the transverse Hall conductivity we find 
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The absorptive part of the longitudinal conductivity 
can be seen in Figs. 0|a) and 0Jb) for the TI and BI 
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FIG. 4. (Color online) Rea xx (Cl) in the (a) TI regime 
(A z = 0.5A SO ) and (b) BI regime (A z = 1.5A so ) for 
B/As O =0.0657T/meV 2 when the chemical potential lies be- 
tween the n — and n = 1 LLs. The first four peaks are spin- 
and valley-polarized while higher features are not. 



regimes, respectively. The transverse Hall conductivity is 
not shown as in the region of spin and valley polarization 
it is simply the negative of the longitudinal conductivity. 
We have used the same parameters as before and included 
a chemical potential of /i — 3.0A SO so that the Fermi level 
lies above all the n = LLs and below all the n — 1 levels. 
This is done so that all the n = LLs are occupied and, 
thus, transitions from n = — 1 to n = are forbidden 
due to the Pauli exclusion principle while all n = to 1 
transitions are permitted. This will occur for any value 
of that is in this region, i.e. |A^-f-| < ±i < Ei(Akx). 
The lower limit, |Axt| 5 is governed by the strength of 
the electric field while the upper bound, Ei(Ak^), can be 
raised by increasing B and/or E z . 

If the chemical potential is situated between n = 
and |n| = 1, we see four absorptive peaks in the longi- 
tudinal conductivity which are spin and valley polarized. 
These peaks result from single spin and valley transi- 
tions as illustrated by the arrows in Fig. [5] While in 
the absence of a magnetic field spin-valley polarization 
has been predicted in the absorptive response to cir- 
cularly polarized ligh t 15 ' 16 , the inclusion of a finite B 
produces four strong, separated peaks of definite spin 
and valley label which can also be observed in the lon- 
gitudinal response. This is to be contrasted with the 
B = response to circularly polarized light where pre- 
dicted valley-spin polarization is not clearly separated 
but contains some admixture of spins and valley o 16 ' 55 . 
Similar to the B = result, spin and valley polariza- 
tion can be seen in the response to circularly polarized 
light for B ^ 0. Figure [5] shows the absorptive part of 
cr±(Cl), where <7±(f2) = <Txx(&) i ia xy (Vt) for right- and 
left-handed polarization, respectively^ 5 -. There remains 
a spin- and valley-polarized quartet of peaks (of double 
weight) in response to right-handed polarization, while 
the quartet is lost in response to left-handed polarization 
and all that remains is the spin- and valley-degenerate 
response at higher energy. In both the longitudinal and 
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FIG. 5. Conductivity for (upper) right-handed and (lower) 
left-handed polarized light. 

circularly polarized responses, the two middle peaks of 
the quartet display a band inversion^ with the transi- 
tion from TI to BI regimes. 

The magnetic field also allows for a tuning of the po- 
sition of the valley-spin polarized peaks over a frequency 
range without having to adjust E z . The onset frequen- 
cies of these peaks is set by the energy difference be- 
tween n = and 1 LLs, namely, fi = £i(A^-^) — Ajf/f, 
e 1 (A K > i ) - Ajf'4., £i(Aif t ) - Aif t and £i(Ajq.) - A^ 
associated with spin up and down at K' and spin up and 
down at K, respectively. Thus, as E\ depends on B, in- 
creasing the magnetic field moves the quartet of peaks 
to higher energy. While the separation between all four 
peaks is not fixed, the separation between the two middle 
peaks is fixed at 2|A^| for A z ^ (all levels are spin 
degenerate when A z = 0) and as a consequence, when 
A z = A so , only the outer peaks of the quartet remain 
spin-valley polarized. This separation is the minimum 
gap between the B = bands and is only controlled by 
E z and the SOC. 

These polarized peaks should be visible in experiment 
as the magnetic and electric field values required to ob- 
serve them are well within experimental limits. Aside 
from B and E z , the determining factor in the onset fre- 
quency of the polarized response is the size of the spin 
orbit gap, A so . For A so = 1.55 meV, the curves shown 
here correspond to a B of ~ 0.15 T and for A so = 7.9, 
B w 4.1 T. In the former case, the splitting of peaks in 
the polarized quartet is ~ 1.5 meV, a resolution easily 
achieved in broadband optics. In this case, the quartet 
sits in the range of 1 — 2 THz where broad band exper- 
iments can be done 5 ^. If B = 1 T, such peaks will shift 
to the far infrared in the range of 17 — 20 meV. If instead 
A so is as large as 8 meV, the splitting of the quartet 
will be about 3-5 meV for B = 1 T and the quartet will 



fall in the range of 13 — 25 meV. Conductivity experi- 
ments on graphene have spanned the range from THz to 
eY 26-29^i with B ^ o T to 18 T^r^^l and resolu- 
tions of order less than a meV; therefore, experiments on 
silicone should observe this polarized behaviour. Exper- 
iments may also provide a measure of A so , for instance, 
by examining when the two middle polarized peaks be- 
come superimposed at the transition between the TI and 
BI regimes (A z = A so ). 

While electron-electron interactions have been a 
source of much discussion in graphen o 58 ' 59 , experimen- 
tal magneto-optics in graphene has been well-described 
by the single-particle picture up to very large magnetic 
fields of ~ 16T— ~ 31 i 59 , hence for modest magnetic fields 
(~1T), electron-electron interactions should not modify 
our main results. In addition, our calculations have been 
based on the assumption of freestanding silicene or sil- 
icene on an insulating substrate where the main effect is 
to provide charge doping, for example, by back gating. 
Currently silicene is fabricated on Ag substrates, but ef- 
forts are underway to find alternative substrates, most 
recently 6 , ZrB2. As in graphene, the choice of insulat- 
ing substrate is not expected to affect qualitatively the 
results presented here. 

Summary: We have predicted the presence of a quar- 
tet of spin- and valley-polarized peaks in the low en- 
ergy electronic density of states which result in spin- and 
valley-polarized absorption lines in the magneto-optical 
longitudinal and transverse Hall dynamical conductivi- 
ties of silicene as well as in the response to circularly 
polarized light. These absorption peaks arise from tran- 
sitions between \n\ = and 1 LLs at both valleys and are 
only present when the chemical potential lies between the 
\n\ = and 1 LLs at the K(K') point. The energy at 
which these peaks occur can be controlled by the strength 
of the magnetic and electric fields. Furthermore, the spin 
and valley labels of the two middle peaks switch as the 
system transitions from a TI to a BI. As the strength of 
the spin- and valley-polarized response is predicted to be 
comparable to what has been observed for graphene and 
to onset at experimentally attainable energies, it should 
be easily visible in experiment. While spin- valley po- 
larization in the absence of a magnetic field has been 
discusse d 15 ' 16 , it can only be isolated with circularly po- 
larized light and may contain some admixture of spin and 
valley. With a finite B, such spin- valley polarization is 
robust, easily separated out and can be observed with- 
out circular polarization. This may prove useful in both 
spintronic and valleytronic technologies as an electronic 
current of definite spin and valley label can potentially 
be isolated. These predictions should hold for other 2D 
crystals with a band gap that can be tuned by an external 
electric field (i.e. germanene). 

We thank J. P. Carbotte for helpful discussions. This 
work has been supported by the Natural Science and En- 
gineering Research Council of Canada. 
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